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ABSTRACT 

We present data from the Supernova Cosmology Project for five high redshift 
Type la supernovae (SNe la) that were obtained using the NICMOS infrared 
camera on the Hubble Space Telescope. We add two SNe from this sample 
to a rest-frame /-band Hubble diagram, doubling the number of high redshift 
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supernovae on this diagram. This /-band Hubble diagram is consistent with a 
flat universe ( ^m, ^^a) = (0.29, 0.71). A homogeneous distribution of large grain 
dust in the intergalactic medium (replenishing dust) is incompatible with the 
data and is excluded at the Sex confidence level, if the SN host galaxy reddening 
is corrected assuming Ry = 1.75. We use both optical and infrared observations 
to compare photometric properties of distant SNe la with those of nearby objects. 
We find generally good agreement with the expected color evolution for all SNe 
except the highest redshift SN in our sample (SN 1997ek aX z = 0.863) which 
shows a peculiar color behavior. We also present spectra obtained from ground 
based telescopes for type identification and determination of redshift. 

Subject headings: supernovae: general - cosmology: observations - cosmological 
parameters 



1. Introduction 



A decade ago, the development of Type la supernova (SNe la) as distance indicators 
began making possible precise measurements of the expansion hist ory, leading to the un- 



derstanding that a dark energy dominates the fate of our Universe (iPerlmutter et al.l Il998 



Garnavich et al. 1998: Schmidt et al. 1998; Riess et al. 1998; Perlmutter et al 



19991). With 



the advent of increasingly larger and more precise SNe la surveys, which aim to elucidate 
the nature of dark energy, i t has beconae clear that systematic uncertainties are a limit- 



Riess et al 



ing factor (Knop et al 



2007 



2003 



Tonrv et al.ll2003l : lAstier et al.ll2006l : IWood-Vasev et al.ll2007b 



Kowalski et al.ll2008l ). In particular, knowledge of SNe intrinsic colors and 



dimming due to extinction by dust in the host galaxy have been identified as the dominant 



system atic effects for the determination of cosmological parameters, see e.g., iNordin et al. 



(120081 ). It is often assumed that dust in distant galaxies resembles the dust in the Milky Way, 
whose properties are well studied and for which maps of dust distribution and composition 
are available. This question, however, remains unresolved. Another source of uncertainty 
is in our knowledge of the SN intrinsic color and its dispersion, which is normally used 



"'^Based on observations made with the NASA/ESA Hubble Space Telescope, obtained from the data 
archive at the Space Telescope Institute. STScI is operated by the association of Universities for Research 
in Astronomy, Inc. under the NASA contract NAS 5-26555. The observations are associated with program 
GO-07850. 
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for determining the amount of reddenings for each supernova. Moreover, the understand- 
ing of the mechanisms leading to the explosion of SNe la or the nature of the progenitor 
systems is still incomplete. Thus, the empirical properties of the standard candles need 
to be continuously compared between distant and nea rby objects. While there is evidence 
for evolution in average SN properties with redshift (IHowell et al.l 120071 ). no intrinsic SN 
evolution has been observed in the spectral and p h otometric properties o f individual object s 



see e.g., 



Nobih et al. (2005): Blondin et a 



Folev et all J2008h : iBronder et all toO^ : 



Riess et al. 



(120071)). 



fl2006h:lGaravini et all (120071 ): Ellis et all (120081 ): 



Dust extinction has a strong wavelength dependence. Milky- Way type dust causes twice 
as much extinction in the i?-band as in the /-band, and a factor of 4 more in the 5-band than 
in the J-band. This suggests that observations in the near-infrared and infrared (A > 0.7 /im) 
are a very powerful tool for minimizing systematic un certainties in super n ova cosmology due 



to the determination of the host galaxy extinction (INobili et al.l l2005l : IWood-Vasey et al. 



2008al ). However, ground-based infrared observations require relatively long exposures and 
are limited by atmospheric emission and absorption. Observations from space have signifi- 
cantly better resolution and less background contamination leading to more precise infrared 
photometry from space than from the ground. 

In this paper, we use Near-Infrared Camera and Multi-Object Spectrometer (NICMOS) 
infrared and Wide Field Planetary Camera 2 (WFPC2) optical data to study the colors of 
five distant SNe la over a broad range of rest-frame wavelengths. In Section [21 we present the 
photometric data. In Section [3], we show the spectra collected for the type identification and 
determination of redshift. We compare the observed optical and infrared colors with those 
of nearby SNe la in Section HI In Section [5l we add two SNe from our sample to the current 
rest-frame J-band Hubble diagram, doubling the number of 2 > 0.1 SNe and extending the 
analysis to z = 0.638. Finally, we discuss our results in Section O 



2. The data set 

We present data of 5 distant Type la supernovae that were obtained with the NICMOS 
mounted on the Hubble Space Telescope (HST). These supernovae, in the redshift range 
0.356 < z < 0.862, were discovered by the Supernova Cosmology Project during ground 



^One normally refers to reddening as an effect of dust extinction, either in the SN host galaxy or in the 
Milky- Way. However, there could be a relation between colors and intrinsic brightness of SNe la that mimics 
an extinction law, though not caused by differential dimming by dust. Thus, we use the term color excess to 
describe the deviation from the average color, irrespective of the underlying mechanism causing the offset. 
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based searches for SNe la conducted with the 4m Blanco telescope at the Cerro Tololo Inter- 
American Observatory (CTIO) in 1997 and 1998. Optical follow up observations were carried 



out f rom the ground and from space and have been presented in previous work flKnop et al. 



20031 ). These consist of multi-epoch observations with the WFPC2 camera on HST in the 
FQ75W and F814W broadband filters and multi-epoch ground-based observations in the 
R and I bands. The HST filter pair and the ground-based filter pair both correspond to 
rest-frame B and V for z < 0.7 and to rest-frame U and B for z > 0.7. 

Infrared images were obtained using the broadband FllOiy filter on the NIC2 camera. 
Given the broad redshift range of the SNe, the FllOVT filter corresponds to different rest- 
frame filters as shown in Figure [TJ 

Some of the NICMOS data were affected by the passage of HST through the South 
Atlantic Anomaly prior to the observations. After visual examination, good images were 
processed w ith a modified versio n of the STScI pipeline. Details of the modifications are 
explained in lFadeyev et al.l (120061 ) and include an improved technique for evaluating per-pixel 
count rate errors and a different method for correcting the cosmic ray (CR) hits. Although 
the standard STScI pipehne also evaluates the count rate errors, the modifications improved 
the statistical correspondence between the errors and fluctuations in the sky level and were 
used to estimate the photometric errors. The standard bias equalization procedure was used 
during the processing. The image s were corrected for the rate dependent non-linearity effect 
jBohhn et al.lboosl : Ide JoneJ [20061 ). 



PSF photometry was done on the processed images with TinyTIM models (IKrist fc Hook 



20041 ). We have used standard star observations to verify that PSF photometry with Tiny- 
TIM is consistent with simple aperture photometry to 1% accuracy. The p rocedure was 
similar to the chi-square fitter used for optical photometry in lKnop et al.l (120031 ). There were 
two main differences due to features of the infrared data, as compared to the CCD-based 
observations. (1) The pixel flux uncertainty was taken from the modified pipeline evaluation. 
(2) Bad pixels, a frequent occurrence in NICMOS arrays, were explicitly omitted from the 
chi-square evaluation. 

For most SNe, the distance to the host galaxy was sufficiently large that the host galaxy 
light could be modeled by a second degree polynomial. In these cases, the reference images 
were used to visually examine the SN location for unusual host galaxy shapes. The only 
exception was SN 1998ay, where an elliptical-like galaxy model was used in a joint fit with 
the final reference images. We have evaluated the systematic uncertainties due to host galaxy 
modeling and CR rejection to be typically at the 1-4% level, although there were a couple 
of cases as large as 10-18%. 
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Because the NICMOS detectors have a count-rate dependen t non-hnearity we correct 



for (often referred to as the Bohhn Effect (IShaw &: de Jong||2008l )). we must apply a corre- 



sponding zeropoint, Znon-iin, to derive Vega magnitudes. Our data were acquired when the 
NICMOS IR arrays were operating at a temperature of 60 K. For this temperature, STScI 
only provides zero points for standard calibrations. The well modeled detector non-linearity 
corrections are based on an operating temperature of 77 K. We assume that the non-linear 
corrections at 60 K are the same as at 77 K§, namely, Zfof - = Zjjf - Zll^_i.^ = 

0.179mag. We obtain Z^^f.^^^ = 22.268. 

A summary of the observations is given in Table [H including the time since i?-band 
maximum, exposure time in seconds, magnitudes in the Vega system, and total uncertainties 
(including both statistical and systematic). Figure [2] shows the infrared image taken at the 
brightest epoch observed. Each image, covering 4x4 arcsecs, is built by co-adding several 
observations for the purpose of displaying the position of the SN with respect to the host 
galaxy. The actual photometry was performed on the individual images. Table [2] reports 
redshifts, 5-band stretch factors. Milky Way reddenings, and color excesses for these SNe. 



3. Spectroscopic confirmation 



Optical spectra of SN 1997ek, SN 1997ez, SN 1998ay, and SN 1998bi were obtained at 
Mauna Kea (Keck II telescope) in December 1997 and April 1998. A spectrum of SN 1998as 
was obtained at the European Southern Observatory 3.6-m telesco pe in March 1998 . All 
supernovae were within 2.5 arcsecs of the host galaxy center. Since iKnop et al.l (120031 ). the 
spectra have been re-analysed and the redshifts re-determined. In general, redshifts are 
measured from host galaxy emission and absorption lines. The one exception is SN 1998bi, 
fo r which the s upern ova was used. Within the uncertainties, the redshifts agree with those 



m 



Knop et al.l (120031 ). 



The extracted spectra, which consist of light fr om both the s upern ova and the host 



galaxy, were fitted with the algorithm described in iHowell et al.l (120051 ). The type, the 



confidence that the supernova is a SN la, and the mean rest-frame spectral age from the five 
best fits are shown in Table [31 The extracted spectra, the host-subtracted spectra, and the 
best matching local SN la are shown in Figure [31 With the exception of SN 1998ay, all of the 
SNe could be identified as either a certain SN la (confidence index of 5) or a highly probable 



•^This is supported by the fact that the NICMOS array response only changed by about 50% after 
installation of the cryocoolcr. On the scale of the non-linearity correction, 0.05 mag per decade of flux 
change, the effect of extra 50% flux nominally results in 0.009 mag correction, a negligible amount. 
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Fig. 1. — Correspondence of the observed R, I and FllOiy passbands (normalized to 1) de- 
reshifted to match the rest-frame optical filters (arbitrarily normalized to 0.6). The dotted 
line show the SNe la spectral template at the time of maximum. 
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Table 1. Summary of Infrared Images 



SN MJD Epoch^ Exp*" Counts^ Vega magnitudes'^ 



SN 1997ek 


50818 


50 


-0 


07 


1026 





483 


± 





055 


± 





014 


23 


058 


± 





126 




50824 


34 


5 


77 


2053 





508 


± 





020 


± 





017 


23 


003 


± 





055 




50846 


73 


28 


16 


1280 





441 


± 





033 


± 





049 


23 


157 


± 





144 




50858 


69 


40 


12 


2560 





253 


± 





017 


± 





017 


23 


760 


± 





101 




50872 


10 


53 


53 


2560 





171 


± 





021 


± 





030 


24 


186 


± 





232 


SN 1997ez 


50818 


64 


5 


34 


2052 





466 


± 





018 


± 





Oil 


23 


097 


± 





050 




50824 


55 


11 


25 


2052 





523 


± 





019 


± 





Oil 


22 


972 


± 





047 




50858 


63 


45 


33 


2560 





231 


± 





017 


± 





007 


23 


859 


± 





087 




50872 


27 


58 


97 


2560 





232 


± 





018 


± 





010 


23 


854 


± 





096 


SN 1998bi 


50922 


51 


10 


70 


3593 





607 


± 





014 


± 





015 


22 


810 


± 





037 


SN 1998as 


50909 


35 


13 


23 


5134 





882 


± 





013 


± 





026 


22 


404 


± 





036 




50921 


27 


25 


15 


5134 





879 


± 





013 


± 





026 


22 


408 


± 





037 


SN 1998ay 


50910 


33 


14 


15 


4105 





554 


± 





017 


± 





012 


22 


909 


± 





040 




50920 


28 


24 


10 


4105 





409 


± 





015 


± 





009 


23 


239 


± 





047 



'^Days since B-band maximum in the observer frame. 
'^Exposure time in seconds. 

'^Counts (ADUs) per second. For c o nversi on to electrons per second one can use the 



Barker fc DahlenI (|2007r ). The first uncertainty shown is statistical 



gain of 5.4 e-/ADU 
and the second one is systematic. 

■^Infrared Photometry. The uncertainty includes both statistical and systematic com- 
ponents, as well as a 2% zeropoint calibration uncertainty. 
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SN 1998 ay 




Fig. 2.— (From top to bottom, left to right) NICMOS images of SN 1997ek, SN 1998as, 
SN 1998bi, SN 1997ez, SN 1998ay and their respective host galaxies. SN 1998ay exploded 
close to the center of its host galaxy. An image of the host of SN 1998ay was taken about one 
year after the explosion, and it has been used for correcting the host galaxy contamination 
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Fig. 3. — Confirmation spectroscopy. With the exception of SN 1998ay, each SN is shown twice. 
The observed spectrum, including light from the host galaxy is plotted in the upper panels. In these 
spectra, the most prominent host galaxy lines, if any, are marked. The red circled crosses mark 
regions of strong telluric absorption. In the lower panels, the binned host galaxy subtracted spectra 
(the black histogram) is shown together with the best fitting local SN la (the blue continuous line). 
The spectrum of SN 1998ay is plotted only once as light from the host galaxy dominates the 
spectrum. 
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one (confidence index of 4) . Tlie spectrum of SN 1998ay is lieavily contaminated with host 
galaxy light so a secure classification based on the spectrum is not possible. However, as we 
show in the next section, the light curve and color of SN 1998ay are consistent with those of 
a SN la. 



4. Supernova colors 



We use the infrared data together with the optical data presented by lKnop et al.l (120031 ) 
to compare photometric properties of distant SNe la to those o f nearby SNe la. Intrinsi c color 
properties of nearby SNe la have been extensively studied by iNobili k. Goobar J2008h . The 
reported color evolution with stretch and supernova epoch are used in the analysis presente d 
here, both as a mod el for the intrinsic colors and for the ii'-corrections (IKim et al.lll996l ). 
Nugent et al.l (120021 ) showed that the latter depends mainly on the supernova color rather 
th an on change s in ind ividual spectral features. Thus we use the spectral ter nplates developed 
bvlHsiao et al.l (120071 ) modified by the color-stretch relation determined in iNobili fc Goobar 
(Eoosl) Defore computing the i^'-corrections. To compare the observed colors of the SNe 
in this data-set with color templates derived for low-z supernovae we need to know three 
parameters: the epochs of our observations with respect to restframe 5-band maximum, the 
stretch s of the rest-frame 5-band lightcurve and the color excess, E{B — V). To obtain 
these parameters, the lightcurves need to be fitted iteratively, since the i^'-corrections needed 
for the fits depend on both s and E{B — V). We use the following procedure: 



Compute i^-corrections using the s = 1 spectral template c orrected for the color excess, 
i.e. warped by the extinct ion law by Cardelli et al.l ( 1989 ) (CCM) with the values of 



E{B — V) for each SN from lKnop et al. I J2003h using a total to selective extinction ratio 
Rv = 1.75, which ha s been found to give a good empirical description of reddening in 
low-z SNla samples ( INobili &: Goobarll2008l . and references therein). 



2. Fit the optical light curves to determine the time of i?-band maximum and the stretch 
factor, s, for each SN. 

3. Compute th e intrinsic colors corres ponding to s following the color-stretch relation 
described in iNobih fc Gooba] j^OOsI) 



Warp the spectral template^ to match the expected colors for the fitted s. We also 
redden the spectral template assuming as a first approximation the color excess E{B — 



^We use a cubic spline interpolation of the ratio between the synthetic photometry and the new photom- 
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V) determined by lKnop et al.l (120031 ) . Thus, we use the modified spectral template to 
compute i^-corrections, while keeping the measured colors unchanged. 

5. Fit the optical light curves again using updated K-corrections 

6. Build the color models using the fitted J-band light curve model interpolated to the 
epochs of the i?-band observations to compute R — I, and to the epochs of the FllOW- 
band observations to compute / — FllOW. 



Note that, in the previous steps we have not corrected the observed photometry for the 
extinction in the Milky Way or in the host galaxy. Instead, we have modified the spectral 
template needed for fitting the data. We use a grid of values of E{B — V) in steps of 0.01 mag 
with a spectral template reddened according to the CCM extinction law. W e assume a value 
of the paramet er Ry = Ay /E {B — V) = 1.75 as in lNobili fc Gooban (120081 ). We note that a 
recent work by iGooban (120081 ) provides a potential explanation for the empirically found low 
average value of the total to selective extinction ratio in Type la supe r novae as bein g a result 
of interaction with dust in the circumstellar environment, (see also IWangl (12005! )). Step 4 
above is repeated until the fit parameters do not change outside the estimated uncertainties. 

Finally, we compute synthetic photometry and build model color curves in the observer- 
frame filters, R — I and / — FllQW , for each SN, reddened by the Milky Way extinction 
( ISchlegel et al.l Il998l ). using the standard value Ry = 3.1. The estimated color excess, 
E{B — V), giving the best fit to the observed optical data, R — I, is reported in Table [21 
These values are then used to redden the color curves / — FllOVT before comparing with 
the observed infrared data. 

Figures HJ- [6] show the comparison between the color curves obtained by synthetic pho- 
tometry and the observed data. The filled symbols show the observed colors, whereas the 
color curves (solid lines) are the expected color evolution. The das hed lines represent the in - 
trinsic dispersion for the given color, determined using the results in lNobili &: Goobarl (120081 ). 
Table [2] summarizes the results of this procedure. We note that the val ues of the stretch 5 
and color excess E{B — V) are slightly different from those reported in iKnop et al.l (120031 ). 
The discrepancies depend mainly on the different intrinsic colors assumed, in particular in 
U — B, and on the spectral template and color-stretch relation used. Since the models as- 
sumed here are based on a larger sample, small discrepancies are not surprising (see also the 



etry at the effective wavelengths. The B-band hght curve bv lGoldhaber et al.l (l200l[) is assumed in o r der to 
compute all the UVRI bands given the colors U - B, B -V, V - R and R- I (see lNobih fc Goobail |2008h 
for details). 
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discussion in lNobili &: Goobarl (l2008l )). We note that most of the differences are within one 
standard deviation from the quoted uncertainties. 

We measured a large color excess, E{B — V) = 0.19 ± 0.04, for the lowest redshift SN 
in our sample (SN 1998as aX z = 0.356). As shown in Figured! the FllOW filter at this 
redshift covers rest-frame R, J, and extends to about 11000 A effectively covering rest-frame 
^-band. Thus, the infrared data explore the red region of this SN spectrum which is less 
affected by reddening. These data offer the possibility to explore intrinsic color properties 
of SNe la despite the larg e color excess in the blue bands. Since the analysis presented 
by (INobili fc Goobarl |2008| ) extends only to the /-band, the colors of the spectral template 
beyond this band are obtained by extrapolation to longer wavelengths. However, the spectral 
templates are quite featureless at these wavelengths and we expect very smooth behavior 
from the J-band to the 2;-band. This gives us confidence about the average color evolution 
within the assumed intrinsic dispersion. We note that the data agree with the expected color 
curves. 

SN 1997ek shows an interesting color behavior. Its optical color is slightly bluer than 
expected and consistent with no extinction, but its infrared colors are substantially redder 
than average, as shown in Figure [61 We note that mis- measurement of the SN flux due to host 
galaxy contamination can be excluded due to the smooth profile and low flux level seen in 
the reference image of the galaxy. Moreover, the time evolution does not follow the behavior 
expected around 20 days after maximum brightness, as shown in the left panel of Figure O 
For this reason the color excess has been estimated using only the filled symbols in the figure. 
Given the different behaviors sh own by this SN, we have compared i ts colors with those of the 
very peculiar nearby SN 2000cx Jhi et al.ll2001al : ICandia et al.lboosi : branch et al1l2004[ l. Pre- 
maximum spectra of SN 2000cx resemble that of the over-luminous SN 1991T and its color 
evolution is quite different from that of a normal Type la. The fact that the pre-maximum 
spectrum of SN 1997ek is well fitted by spectroscopically SN 1991T-like SN 1999aa, together 
with its peculiar colors, has prompted us to attempt a comparison of these two SNe. We 
have warped the spectral template with UBVRI photometry of SN 2000cx before computing 
synthetic photometry with the passband filters used for the observations of SN 1997ek. The 
resulting color curves are shown as dotted lines in Figure O The agreement is improved in 
the infrared around maximum, but not at later times where SN 2000cx has colors similar 
to those of a normal SN la. Moreover, the dotted curve does not fit the observed data in 
the optical. We note that, if only the optical data at maximum were available, one would 
conclude that this SN is quite normal. The peculiar color behavior of SN 1997ek could be 
an indication of an intrinsic peculiarity of this SN. We note that the confidence index for the 
spectral typing of SN 1997ek is 4. Although it is not as secure as candidate that has clear 
SNe la features, such as Si II or S II, it is unlikely to be a SN of a different type. 
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The value of E{B — V) determined from the optical data gives a good agreement between 
the model and observed infrared color, giving a value of the ~ dof for all cases except 
SN 1997ek. We take the good agreement found in the case of SN 1998ay as a confirmation 
that this is indeed a Type la SN. 

We compare the total distribution of the residuals of our sample with the expected 
distribution of nearby SNe. We normalize the residual of each data point from the expected 
color by dividing by the total uncertainty. The total uncertainty is calculated as a quadratic 
sum of the measurement uncertainty and the intrinsic dispersion: 

Figure [7] shows the pull distribution for all SNe in each observed color. The shaded area in- 
clude all SNe except the peculiar SN 1997ek. This is consistent with a Gaussian distribution, 
with zero mean and a = 1. Since the intrinsic dispersion used is based on measurements 
of nearby SNe, we can conclude that the same average color dispersion is observed in the 
distant sample. 



4.1. SN 1999ff and SN 2000fr 



Nobili et al.l (120051 ) presented ground-based infrared observations for one SN la, SN 2000 fr 
at redshift z = 0.543. This SN together with SN 1999fr at ^ = 0.455 JTonrv et al. l l2003l ) 
was compared to nearby SNe la and used to build a rest-frame /-band Hubble diagram. We 
present the result of a re-analysis of these data, using the same technique applied to our 5 
SNe in the previous section. Figure [8] shows optical and infrared color evolutions for both 
SNe. Table H] summarizes the results. 

We note that the color excess giving the best fit to the optical data is consistent with 
no reddening for both SNe. The infrared data also agree with the expected behavior within 
the uncertainties. 



5. /-band Hubble diagram 

Two SNe in our sample are at a redshift where the infrared observations correspond 
to the rest-frame /-band (see Figured]). These are SN 1998as {z = 0.356) and SN 1998ay 
[z = 0.638). We note that SN 1998as has the largest extinction in our sample and SN 1998ay 
does not have a confirming spectrum, but is fully consistent with a SN la due to its colors 
and light curve shapes. We used these SNe to extend the /-band Hubble diagram built by 
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Fig. 4. — Observer-frame color for each SN. The color curves are obtained by synthetic 
photometry of the spectral templates (see text for details). The dashe d lines show the 
intrin sic dispersion around the color curves measured on nearby SNe la (INobili fc Goobar 
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Fig. 5. — Observer-frame color for each SN. The color curves are obtained by synthetic 
photometry of the spectral templates (see text for details). The dashe d lines show the 
intrin sic dispersion around the color curves measured on nearby SNe la (INobili fc Goobar 

mi). 
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Fig. 6. — Observer-frame color for SN 1997ek (z=0.863). The color curves are obtained 
by synthetic photometry of the spectral templates warped to match the intrinsic colors of 
SNc la (solid line) and warped to match the colors of the peculiar SN 2000cx (dotted line). 
Sec text for details. The open symbol has not been used for computing the color excess. 
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Table 2. Supernova Name, Redshift, -B-band Stretch Factor, Milky Way Reddening, and 

Color Excess. 



SN 



SB 



E{B - V) 



MW 



E{B - V) 



SN 1997ekb 0.862 

SN 1997ez 0.776 

SN 1998as 0.356 

SN 1998ay 0.638 

SN 1998bi 0.75 



1.059 ± 0.071 0.042 

1.094 ± 0.038 0.026 

0.931 ± 0.018 0.037 

1.050 ± 0.045 0.035 

0.975 ± 0.031 0.026 



-0.04 ±0.05 
0.20 ±0.08 
0.19 ±0.04 
0.09 ±0.09 
0.16 ±0.04 



'^Milky Way reddening given by ISchlegel et al.l (119981 ) 

^'The color excess for SN 1997ek is estimated excluding the data 
point at day ~ 20. See text for details. 




Fig. 7. — Histograms of the residuals of all data from the expected color curves (solid line) 
and for all SNe but SN 1997ek (shaded area). For the latter we measure a = 0.95, for N = 25 
(left panel) in the optical; and a = 0.97, for iV = 9 (right panel) in the infrared. See text 
for details. 



Table 3. Summary of Spectroscopic Observations. 



SN 


Telescope 


Instrument /Setup* 


Date(UT) 


Exposure'^ 


Best Match 


Host Lines 


Epoch'^ 


Index 


Classification 


SN 1997ek 


Keck II 


LRIS/400/8500/1 


02/01/1998 


3x1500 


SN 1999aa -3 days 


[0 II],H,K 


-1.8 (4.1) 


4 


SN la 


SN 1997ez 


Keck II 


LRIS/300/5000/1 


31/12/1997 


3x1800 


SN 1994D +2 days 


[0 II] 


+0.4 (6.4) 


4 


SN la 


SN 1998as 


ESO 3.6m 


EFOSCII/R300/1.5 


28/03/1998 


2x1800 


SN 1989B -1 day 


H/3, [0 III], Ha 


-1.8 (4.3) 


5 


SN la 


SN 1998ay 


Keck II 


LRIS/400/8500/1 


03/04/1998 


2x1200 


No match 


[0 II], H/3, [OIII] 




2 


SN 


SN 1998bi 


Keck II 


LRIS/400/8500/1 


01/04/1998 


1x1500 


SN 1994D -4 days 




-3.4 (3.9) 


5 


SN la 



Tor the LRIS observations, 400/8500/1 means that the observations were done with the 400/8500 grating and the 1" slit. For the EFOSC II observations, 
R300/1.5 means that the observations were done with the R300 grism and the 1'.'5 slit. 

^Exposure time expressed in seconds 

'^Average rest-frame epoch of the five best fits, with the rms value reported in parenthesis I 

I— ' 
00 

I 
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Fig. 8 — Observed-fram e color for SN 1999ff flTonrv et al.ll2nn3l ) and SN 2000fr flKnop et al. 
20031 : iNobili et al.ll2005l ). The color curves are obtained by synthetic photometry of the 
spectral templates (see text for details). The dash ed lines show the intr insic dispersion 
around the color curves measured on nearby SNe la (INobili &: Goobarll2008l ). 
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f lNobili et all l2005[ ) which originally included 2 z > 0.1 SNe: SN 2000fr at z 
SN 1999ff at 2 = 0.455. 



0.543 and 



We computed /T-c orrections for all four SN e from the observed infrared filter to rest- 
frame /-band follow ing iNobili fc Goobarl (120081 ). Using the J-band templates of 42 nearby 
SNe la developed in lNobili et al.l (120051 ) we fit the peak brightness of the high-redshift SNe. 
Figure [9] shows the template fits that satisfied < Xmm + where Xmm is the of the best 
fitting template. The fitted Imax giving the minimum is reported in Table O However, 
since several templates correspond to a similar value of the x^? we estimated the uncertainty, 
A, as the rms of the peak magnitudes that satisfied x^ < Xmm + 1- 

Note that, given the different way of computing i^-corrections, th e results obtained f or 
SN 2000fr and SN 1999ff are slightly different than what was found in lNobih et all (|2005[ l. 



The peak magnitudes have been corrected for the width-luminosity relationship in the 
/-band and the Milky Way absorption, A^"^ . The host galaxy extinction correction, Aj"''*, 
was estimated from the measured color excess and assuming the CCM extinction law with 
Ry = 1.75: 



mf = mj + aj{sB - I) - A^"'' - A''''^ (1) 
The effective m agnitudes of the fou r high redshift SNe together with the effective magnitudes 



of nearby SNe (iNobili et al.l l2005l ) have been used to build the rest-frame /-band Hubble 
diagram shown in Figure [TOl A value of a/ = 1.18 ± 0.19 was fitted for the low-redshift 
sample and used for all SNe. The uncertainties reported in Table [5] are represented by the 
inner error bars, and the larger error bars are obtained by adding an intrinsic dispersion 

Table 4. Parameters of the Additional SNe, Redshift, /?-band Stretch Factor, Milky Way 

Reddening, and Color Excess. 



SN^ 


z 


SB 


E{B — V)mw^ 


E{B - V) 


SN 1999ff 


0.455 


0.82 ± 0.05 


0.025 


0.06 ±0.11 


SN 2000fr 


0.543 


1.054 ± 0.012 


0.030 


0.01 ±0.04 



'^From lNobih et all (120051 ) 
^Given by ISchleeel et all Jl998h 
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Fig. 9.— Best fits to rest-frame /-band observations for SN 1998as, SN 1998ay, SN 1999ff 
and SN 2000fr. All templates giving a < xLin + 1 shown together with the best fit 
template (thicker red solid line). The best fit templates are listed in Table [5] 
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equal to 0.11 mag as measured on the nearby SNe. The rms dispersion of the high-redshift 
SNe is 0.15 mag. 

The sohd hne represents the best fit to the rest-frame /-ban d nearby data for the 
concordance model with fixed VLm = 0.29 and = 0.71, as foun d bvlKowalski et ahl (120081) 



in the S-band. The single fitted parameter, A^/, is defined as in iPerlmutter et al.l (119971 ): 



A4/ = M/-51ogi/o + 25 



(2) 



where M/ is the /-band absolute magnitude for a S-band stretch = 1 supernova. The 
value fitted is Adi = —3.19 ± 0.02. The dotted line in Figure [TOl represents a cosmological 
model without dark energy {Qm = 1 and Q\ = 0) while the dashed-dotted line represents 
a cosmological model without dark matter [Qt = 0.04 and Qa = 0.96), where the matter 
density is equal to the baryon matter density only. Also plotted is the model of a universe 
without dark energy in the presence of a homogeneou s distr ibution of large grain dust in 
the intergalactic medium (dashed line). iGoobar et al.l (120021 ) hypothesized an ad-hoc dis- 
tribution of dust that is able to explain the dimming observed in the SNe la /?-band peak 
magnitudes. According to this model, the comoving density of the dust increases with cosmic 
time until z=0.5. This "replenishing dust" is created at the same rate as universal expansion. 
By construction it is very diffic ult to constrain t his kind of dust model using the /?-band 
Hubble diagram alone (see e.g., iRiess et al.l 120071 ). We are able to exclude this model with 
high confidence with the rest-frame /-band Hubble diagram, as shown in Table El All values 
of computed with the four high-z SNe for each model tested are also reported in the table. 

To see the effect of different values of Ry on the Hubble diagram, we repeated the 
analysis in the case of a standard Milky Way type dust also in the host galaxy, with Ry = 3.1. 
The results are summarized in Table El It is worth noticing that the constraint put on the 
replenishing dust model are not equally stringent in this case. The residuals of the high 
redshift SNe are shown in the inset in Figure [TOl We note that the dispersion of the data 
set is smaller in this particular case for Ry = 3.1 than for Ry = 1.75. 



We repeat the cautionary note in lNobili et al.l (120051 ). The different method applied to 
determine the peak magnitudes for the nearby and distant SNe la might introduce systematic 
uncertainties which are difficult to quantify. Moreover, many (but not all) of the templates 
used for fitting the high-redshift SN are built on the nearby supernovae in the Hubble diagram 
leading to a possible source of correlation. 
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Table 5. Results of the Rest-Frame /-band Light Curve Fits for Two of the SNe in our 
Sample, Together with SN 2000fr and SN 1999ff. 



SN 


z 


^max 




Best Template 


Templates'^ 


SN 1998as 


0.356 


23.004 


0.146 


1998ab 


1993H,1994D,1997E,1998V,1999cl 


SN 1998ay 


0.638 


24.175 


0.048 


1994T 


1992bp,1997dg 


SN 2000fr 


0.543 


23.495 


0.033 


1992bc 


1992bh,1994ae,1996bl,1999gp 


SN 1999ff 


0.455 


23.359 


0.058 


1992bc 


1992bh,1999gp,1996bl,1999aa, 



1996C,1995D,1992bh 



^Peak magnitude obtained by the best fit template. 

'^The rms magnitude of the best-fit templates that satisfied < Xmm + 1- See text for 
details. 

^Best fit templates giving a < Xmin + 1- ^^e text for details. 



Table 6. with 4 d.o.f. for the Cosmological Models in Figure [TO] 





^2a 


p{x^ < xir 


^2b 




(0.29,0.71) 


7.04 


0.87 


0.78 


0.06 


(0.04,0.96) 


6,04 


0.80 


6.21 


0.82 


(1,0) 


54.06 


1.00 


25.71 


1.00 


(1, ^^dust 


23.64 


1.00 


7.75 


0.90 



'^Results obtained assuming Ry = 1.75 
Results obtained assuming Ry = 3.1 



-24- 




Fig. 1 0. — Rest-frame /-band Hubble diagram. The open symbols are data from lNobili et al.l 
( 120051 ) ■ The solid symbols are 2 high-z SNe from this work. The solid line is the best fit 
to the nearby data for the concordance model with fixed Qj^j = 0.29 and Q\ = 0.71. The 
dotted line is a fiat universe model with no cosmological constant, while the dashed line is 
the same model in presence of large grain dust in the intergalactic medium. The dashed- 
dotted line is a fiat universe with no dark matter, Qm = 0.04, and = 0.96. The inner 
error bars represent the uncertainties reported in Table [5l while the outer error bars are 
obtained adding an intrinsic dispersion equal to 0.11 mag. The inset shows the effect of 
using Ry = 3.1 for host galaxy dust extinction. 
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6. Summary and discussion 



This data set is unique in terms of wavelengtli coverage and redshift. In fact, only a 



handf ul of SNe la at these redshifts have multi-epoch observations in the near-IR. iRiess et al. 



(120071 ) presented NICMOS observations of SNe la at z > 1. However, the majority of the 
data correspond to rest-frame B and V and are mostly at later epochs between 5 and 20 
days after maximum brightness. 

In Section[3l we presented the spectra that were used for SN typing and for measuring the 
redshift. All SNe were securely identified as Type la except SN 1998ay which has a spectrum 
suffering from severe host galaxy contamination. SN 1998ay shows all other characteristics 
of a Type la SN. The observed colors are in perfect agreement with the models in the optical 
and in the infrared, providing strong evidence that this is a Type la SN. 

The high quality of the infrared data allowed a detailed comparison of the observed 
optical-IR colors of these high-redshift SNe la to the colors of nearby SNe la. This was 
presented in Section |H Optical colors were used to determine the color excess in rest-frame 
E{B — V). The color excess correction to be applied to the infrared was determined assuming 
the COM extinction law. We note that the small amount of absorption in the infrared 
compared to the optical bands and the relatively large uncertainties in the measurements 
undermine the possibility of determining Ry- This is even more severe for those SNe with 
small color excesses suc h as SN 1998ay. For this reason we have assumed a nominal standard 
value of Ry = 1.75 as in lNobili &: Goobarl (120081 ). Assuming different values of Ry, including 
the commonly assumed value of 3.1, lead to estimates with similar color-excess. We note a 
generally good agreement with the expected color evolution in both R — I and in I — FllOW . 

SN 1997ek is the most distant SN in our sample and the only SN with unusual colors. It 
is found to have a normal optical color around maximum light corresponding to rest-frame 
U — B. However, it is significantly redder than average in the observer-frame I — FllOVT, 
corresponding to rest-frame B — VRel, where VRel is a combination of V, R, and a small 
fraction of /. This behavior is not compatible with a reddening process which by definition 
would affect shorter wavelengths more than longer ones. Thus, we take it as an indication 
of a peculiarity in the intrinsic color of this SN. We fo und no counterpart of this object in 
the nearby sample analyzed by lNobili fc Goobarl (120081 ). 



The spectrum of SN 1997ek at 3 days pre-maximum is similar to that of the peculiar 



SN 1 999aa, which is spectroscopically similar to the over-luminous SN 1991T (IGaravini et al. 



2004f ). Including all 1991T-like, 1991bg- like and 1986G- like SNe, about 36% of all SNe la in 
the nearby sample are peculiar SNe la (ILi et al.ll2001bl). However, only a small fraction of 



peculiar SNe la are discovered at high redshift (see e.g.. iMatheson et al.ll2005l : iBlondin et al. 



- 26 - 



20061 : IGaravini et al.l 120071 : iBronder et al.ll2008l ). On t he other hand who lly unknown SN 
subtypes are discove red at high-z, e.g., SNLS-03D3bb (IHowell et al.l l2006l ) and SCP 06F6 
( iBarbary et al.ll2009l ). The different rate of peculiar objects discovered could be taken as 
an indication of evolution of SN properties with redshift. However, other explanations are 
also possible. For instance, magnitude-limited surveys could fail to detect sub-luminous 
SNe or over-luminous SNe suffering high dust extinction. The lower signal-to-noise ratio in 
the observed spectra of high-z SNe la also lirn its the possibility of discovering spectroscopic 
peculiarities (see discussion in iLi et al.ll2001bl ). Multi band photometry offers a very useful 
tool for discriminating between SN subtypes by allowing a comparison of colors which can 
be measured with higher signal-to-noise ratios. Identifying properties that distant SNe have 
in common with nearby SNe, even if peculiar, supports the validity of SNe la as distance 
indicators to probe the composition of the universe. Future large SNe la datasets with 
multi-band light curves (such as the SNLS) will help to determine whether other SNe like 
SN 1997ek exist at high-z. 

Another intriguing possibility is that this peculiarity is an effect of varying metallicity. 
This would have a small effect on the optical peak bright ness of the supernova but would 



signifi c antly impa c t the UV spectral energy distribution (iHoefiich et al.l Il998l : iLentz et al. 



2000). lElhs et al.l fl2008f ) measured the UV fiux in 36 high redshift SN e in th e SNLS sam- 
ple and did not find significant evolution. More recently, ISauer et al.l (120081 ) showed that 
enhanced Fe and Ti/Cr abundances result in an increase of fiux in the UV. In particular, 
an increase of iron abundance affects the region bluer than 3000 A and leads to a variation 
in the [/-band magnitude of up to 0.5-0.7 mag . We n ote that the F336W pass-band filter 
analyzed in Figures 4 and 5 of the ISauer et al.l (120081 ) analysis match our observed i?-band 
for SN 1997ek. We tested the hypothesis of increased fiux in the [/-band due to metallicity 
for SN 1997ek with Monte Carlo simulations. The spectral template was first modified by 
an increment in fiux corresponding to a decrement in the blue end of the spectrum, AU. 
Color excess corrections to the spectrum were then applied for a range of E{B — V) val- 
ues with Rv = 1.75 before computing the synthetic photometry. Both optical and infrared 
data were used to find the best match to the synthetic color curves: mag 
and E{B — V) = 0.5~^q'1 mag with = 6.26 for 8 dof. The model clearly provides a very 
good fit to this pa rticular supernova. Thus, we have utilized one of the models developed by 
Sauer et al.l (120081 ) . for SN 2001ep, giving a similar AU, to modify the colors of SN 1997ek. 
The model is obtained by increasing the iron abundance by +2 dex. We obtained the best 
fit to our data, for E{B ~V) = 0.40 ± 0.05 mag, for = 13.8 for 9 dof. We note that such 
value for the color excess, would imply a large correction to the 5-band peak magnitude, 
that would in turn move this SN about 1 mag away from the best cosmology fit in the Hubble 
diagram. 
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Finally, we used the infrared observations of SN 1998ay and SN 1998as to build a rest- 



frame /-band Hubble diagram. We have added these two SNe to the sample from lNobili et al. 



(120051 ) to extend the Hubble diagram up to ^ = 0.638, doubling the high redshift sample. All 
four high redshift SNe were treated consistently to fit the /-band peak magnitudes. We found 
the data to be consistent with a A dominated fiat universe, with {flu, ^a) = (0.29, 0.71), for 
both values of Ry assumed. The smaller dispersion in the /-band and the smaller extinction 
correction uncertainties allow us to exclude the presence of a homogeneous distribution of 
large grain dust in the intergalactic medium at the five standard deviation level with only 4 
SNe la for the case of Ry = 1.75. This conclusion is severely weakened for Ry = 3.1. We 
note, however, that the data set used in this paper is very inhomogeneous. Moreover, we 
applied a different method for fitting the light curve peak magnitude for the nearby and the 
distant SNe which could be a source of systematic uncertainties. 

Multiband observations such as those provided by the proposed SNAP satellite will be a 
powerful method for sub-typing of SNe la and for determining precise distances. Observations 
of hi gh redshift SNe la i n the near infrared will extend the rest-frame /-band Hubble diagram 



(e.g.. 



Burns et al.ll2006l ). Already proven to be an excellent test for systematic uncertainties, 
this approach will continue to compliment the more extensively used rest-frame SN B-hand 
Hubble diagram. 
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